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The linear and nonlinear transport properties of the single-electron transistor at the degeneracy 
point are investigated for the case of weak single-mode tunnel junctions. Two opposing scenarios 
are considered, distinguished by whether or not electrons can propagate coherently between the 
two tunnel junctions. Each of these two scenarios corresponds to the realization of a different 
multichannel Kondo effect — the two-channel Kondo effect in the case where coherent propagation 
is allowed, and the four-channel Kondo effect in the absence of coherent propagation. A detailed 
analysis of the linear and nonlinear conductance is presented for each of these scenarios, within a 
generalized noncrossing approximation for the nonequilibrium multichannel Kondo Hamiltonian. A 
zero-bias anomaly is shown to develop with decreasing temperature, characterized by the anomalous 
power laws of the multichannel Kondo effect. A scaling function of the differential conductance with 
V/T {V being the applied voltage bias, T the temperature) is proposed as a distinctive experimental 
signature for each of these two scenarios. In the absence of coherent propagation between the 
leads, and for asymmetric couplings to the two leads, a crossover from four-channel to two-channel 
behavior is manifested in a vanishing zero-temperature conductance, and in a nonmonotonic voltage 
dependence of the differential conductance for small asymmetries. 



PACS numbers: 73.40. Gk, 73.23.Hk, 72.f5.Qm 
I. INTRODUCTION 

The single-electron transistor — a small metallic or 
semiconducting quantum box connected to two separate 
leads — is among the basic elements of mesoscopic de- 
vices. Due to the finite energy barrier for charging the 
box with a single electron, the so-called Coulomb block- 
ade, ^'^ charge inside the box is nearly quantized at low 
temperatures for weak tunneling between the box and the 
leads. As a result, transport through the box is strongly 
suppressed unless the two lowest lying charge configura- 
tions in the box are tuned to be degenerate. As a function 
of gate voltage, the conductance thus shows a sequence 
of narrow peaks, each corresponding to the crossing of 
the ground state from n to n -t- 1 excess electrons inside 
the box. 

At the degeneracy points, the system is subject to 
strong charge fluctuations. The corresponding low- 
energy physics is governed by the non-Fermi-liquid fixed 
point of the multichannel Kondo effect, where the two 
degenerate charge configurations in the box play the role 
of the impurity spin. Which multichannel Kondo effect 
is realized depends on microscopic details such as the 
number of transverse modes in the junctions, and the na- 
ture of electron transport inside the box. Two opposing 
scenarios were considered to date for the single-electron 
transistor, distinguished by whether or not electrons can 
propagate coherently between the two tunnel junctions. 
Focusing on wide tunnel junctions, Konig et al.^ consid- 
ered the case where electrons can propagate coherently 
between the two leads (see also Refs. 7 and 8). Extending 
the work Grabert^ to nonequilibrium transport, these au- 
thors analyzed in detail all second-order contributions to 
the current in the dimensionless tunneling conductance. 



obtaining good agreement with experiment. However, 
based on perturbation theory, this approach breaks down 
near the degeneracy points, where transport is governed 
at low temperature by the strong electronic correlations 
of the multichannel Kondo effect. 

An alternative scenario was considered by Furusaki 
and Matveev,^^ and later by Zarand et al}^ and Le Hur 
and Seelig.^'^ Noting that elastic cotunneling is strongly 
suppressed at temperatures above the level spacing,^'' 
these authors omitted altogether coherent electron trans- 
port between the leads, by coupling each lead to in- 
dependent conduction-electron modes within the box. 
For symmetric single-mode junctions, the resulting low- 
temperature physics is governed at the degeneracy point 
by the four-channel Kondo effect, "'^^ in contrast to the 
two-channel Kondo effect that takes place when elec- 
trons can propagate coherently between the leads. Any 
asymmetry in the coupling to the two leads drives the 
system away from the four-channel fixed point to a two- 
channel fixed point, where one lead is decoupled from 
the box. Consequently, the zero-temperature conduc- 
tance vanishes, as shown by Furusaki and Matveev in 
the limit of both a large asymmetry and strong tunnel- 
ing to one lead^^ (i-e., a nearly open tunneling mode). A 
quantitative description of the temperature and asymme- 
try dependence of the conductance in this case remains 
lacking. 

Despite considerable efforts, the understand- 
ing of the low-temperature transport at the degeneracy 
points is far from complete for either scenario. For weak 
single-mode tunnel junctions, there is no quantitative 
theory for the temperature dependence of the conduc- 
tance in the Kondo regime, while the nonequilibrium 
differential conductance is practically unexplored in this 
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FIG. 1. Schematic description of the physical system. A 
quantum box is coupled capacitively to a gate, and con- 
nected by weak tunneling [tL and tR) to two metallic leads: 
a left lead and right lead. The charge inside the box is con- 
trolled by varying the gate voltage Vg, which fixes the electro- 
static potential within the box. A drain-source voltage bias 
V = Vr — Vl is applied across the transistor, which drives a 
current through the device. 



regime for either scenario. The goal of this paper is 
to provide a detailed analysis of the linear and nonlin- 
ear transport at resonance, for weak single-mode tunnel 
jimctions. Both scenarios where electrons either can or 
cannot propagate coherently between the two junctions 
are considered. Our aim is to provide a host of signa- 
tures that can be used to experimentally discern the two 
pictures, and to detect which multichannel Kondo effect 
is realized in actual systems. "'^^ 

To this end, we employ the noncrossing approxima- 
tion^^ (NCA) . The NCA is a self-consistent perturbation 
theory about the atomic limit. Originally designed to 
study dilute magnetic alloys, this approach was success- 
fully applied to the out-of-equilibrium Kondo effect both 
for the single-channel-*^^"^" and two- channel^ -'^ Anderson 
impurity model. Recently, the NCA was generalized to 
the multichannel Kondo spin Hamiltonian with arbitrary 
spin-exchange and potential-scattering couplings. Simi- 
lar to the NCA formulation of the multichannel Anderson 
model,^^ the Kondo-NCA (KNCA) correctly describes 
the low-energy physics of the multichannel Kondo model, 
reproducing the exact non-Fermi-liquid power laws and 
logarithms of the multichannel Kondo effect. This should 
be contrasted with the single-channel case, where the 
NCA (and KNCA) fails to describe the Fermi-liquid fixed 
point. -'^^ Here we extend the KNCA approach to the 
nonequilibrium case, and apply it to the single-electron 
transistor. 

Using the KNCA, we find that a zero-bias anomaly de- 
velops in the current characteristics of the single-electron 
transistor, featuring the anomalous power laws of the cor- 
responding multichannel Kondo effect. The latter power 
laws are manifest in the temperature and voltage depen- 
dences of the low-temperature, low-bias differential con- 



ductance. Similar to the case of two-channel Kondo scat- 
tering off nonmagnetic two-level timneling systems, ^^'^^ 
these power laws are best revealed in a scaling plot of the 
differential conductance versus eV/ksT, with V being 
the applied voltage bias. Such scaling plots are proposed 
as a distinct experimental signature of the appropriate 
multichannel Kondo effect. In the absence of electron 
propagation between the two leads, a crossover from four- 
channel to two-channel behavior is found in the differen- 
tial conductance for asymmetric couplings, in accordance 
with the picture of Furusaki and Matveev.^^ 

The remainder of the paper is organized as follows. In 
sec. II we present the single-electron transistor, and dis- 
cuss the two alternative models under consideration. The 
nonequilibrium KNCA is then formulated in sec. Ill, fol- 
lowed by the derivation of the current in sec. IV. Our re- 
sults for the linear and nonlinear transport are presented 
in sec. V, and discuss in sec. VI. 



II. THE SINGLE-ELECTRON TRANSISTOR 

A. Basic model 

The physical system under consideration is shown 
schematically in Fig. 1. A metallic island, or quantum 
box, is connected by narrow point contacts to two sepa- 
rate leads, a left lead (L) and a right lead (R). A drain- 
source voltage bias is applied across the device, which 
sets a chemical-potential difference between the leads: 
jJ-L — IJ-R = eV. Here — e is the electron charge. Model- 
ing each lead by N noninteracting one-dimensional con- 
duction modes, the Hamiltonian of the two independent 
leads reads 



N 



a=L,Rn=l k,cr 



(1) 



where c 



t 

ockna 



creates a conduction electron with wave 



number k and spin projection cr in the nth mode of lead 

a, and e^kn are the corresponding single-particle ener- 
gies, measured relative to the chemical potential on that 
lead. 

For the quantum box, one has to consider also its 
charging energy. Setting the Fermi energy of the box as 
our reference energy for the single-particle levels inside 
the box, the excess number of electrons on the metallic 
island is described by the operator 
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m=l q,<T 



(2) 



Here, similar to our notation for the leads, c^Bqma 
ates an electron with wave number q and spin projection 
a in the mth mode of the box (M independent conduc- 
tance modes are taken within the box), and eBqm ai"© the 
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corresponding single-particle levels. The latter levels are 
assumed to be sufficiently dense such that a continuum- 
limit description can be used. The Hamiltonian of the 
isolated box is thus given by 

M 2 
Wfeox = ^2 5Z ^Bqm.C^Bqma'^Bqma + Ec — Nc^ , 

m=l q,(7 

(3) 

where Eq is the charging energy of the box. and Nq is 
the classical number of excess electrons inside the box. 
These two quantities are related to the capacitances of 
the gate and of the left and right junctions through Ec = 
e'/2{Cg + CL + CR) and Nc - [CgVg + CLVL+CRVR) / e, 
where Vg, Vl, and Vr are the voltages applied to the gate 
and to the left and right leads, respectively (see Fig. 1). 

The quantum box is coupled to the leads by weak tun- 
neling, described by the tunneling Hamiltonian 

Utunnel = ^ ^nm {c^akna^Bqma + H.C.j . (4) 

a=L,R. n.m k,q,a 

For simplicity, both the k and q dependences of the tun- 
neling matrix elements and have been neglected 
in Eq. (4). The full Hamiltonian of the system reads 

^ — "Pleads "i" "^box "^tunnel * 

B. Mapping onto the multichannel Kondo problem 

Our objective is a detailed quantitative description of 
the linear and nonlinear transport at resonance, when 
two neighboring charge configurations are degenerate 
within the box. To this end, let us focus on the vicinity 
of a particular degeneracy point, Nc = n + ^ with n an 

integer, separating the Nb = n and Nb = n + 1 charge 
configurations. 

As emphasized by Matveev,^^^ the low-tcmpcrature, 
low-bias physics is governed at resonance by the 
intermediate-coupling fixed point of the multichannel 
Kondo effect. To formulate the connection between the 
two problems, we proceed along the lines laid out by 
Matveev.^ Labeling the deviation from the degeneracy 
point by ^A'' = Nc — n — ^, we concentrate on \5N\ <C 1 
and e\V\,kBT <C e^/Cs, such that all charge configura- 
tions in the box other than Nb = n and Nb =n+l are 
energetically inaccessible. One can formally remove all 
higher energy charge configurations in the box by means 
of two projection operators, Pn and Pn+i, which project 
onto the Nb = n and Nb = n+1 subspaces, respectively. 
Alternatively, one can carry out the projection onto the 
low-energy subspacc by introducing a spin-^ isospin op- 
erator S, and identifying 

Pn+l — Pn < > 252, (5) 

-Pn^Bqmcr-Pn+1 ■* ^ ^Bqmcr^ ) (6) 
Pn+lC^Bqma^n •* > ^^Bqma^^ ■ C^) 



Omitting an E'c'-dependent reference energy, the result- 
ing low-energy Hamiltonian reads 

Weff = ^ ^ {iakn + ^a)Cakncr^o:kna 

a—L^R k,n.(T 

+ ^BqmC^Bqma^Bqma - 2EcSNSz 

q,m,a 

+ E E {inm<ikna(^BqmaS- +-R.C.} . (8) 

cx=L,R k,q,n,Tn,a 

Strictly speaking, Eqs. (5)-(8) are subject to the con- 
straint Nb — Sz = n + 1/2, refiecting the fact that the 
isospin 5 and the conduction-electron operators c^^^^ 
are not independent degrees of freedom. However, since 
the isospin dynamics in the Hamiltonian of Eq. (8) is not 
sensitive to the precise number of conduction electrons 
inside the quantum box, this constraint can be conve- 
niently relaxed. We therefore regard hereafter the isospin 
S and the c^Bqma operators as independent. 

The Hamiltonian of Eq. (8) has the form of a plan- 
ner multichannel Kondo Hamiltonian. Here h = 2EcSN 
plays the role of a magnetic field, while the equivalence 
of the two spin orientations guarantees the existence of 
at least two identical conduction-electron channels. The 
total number of independent conduction-electron chan- 
nels depends, however, on the microscopic details of the 
tunneling matrix elements. In this paper we consider 
two particular scenarios, corresponding to the two- and 
four-channel planner Kondo Hamiltonian. 

C. Coherent propagation versus no coherent 
propagation between the two junctions 

For sufficiently narrow point contacts, only a single 
mode weakly couples each lead to the quantum box. Fo- 
cusing on this case, we consider two different scenarios: 
one by which both leads couple to the same single mode 
within the box, and the other whereby each lead is cou- 
pled to a different mode within the box. 

The former scenario is just a single- mode version of the 
model of Konig et al.f' whereby electrons can propagate 
coherently between the two leads. Its corresponding low- 
energy effective Hamiltonian is obtained by setting N = 
M = 1 in Eqs. (1) (4), thus omitting the mode indices n 
and m. The resiilting Hamiltonian takes the form 

n..= E E( 

£afc + Ma)c]j/jo-'^afe<^ 5^ ^BqC'Bqa^Bqcr 

a=L,R k,a q,<T 

+ E Y.^'-i'^'^ccka'^BqaS- +ii.c]-hSz, (9) 
(x=L,R k,q,CF 

where II and Ir (both taken to be real) 

are the tunneling matrix elements for the left and right 
junctions, respectively. 
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The second scenario is a straightforward adaptation 
to weak single-mode tunnel junctions of the model in- 
troduced by Furusaki and Matveev,^^ and later used by 
Zarand et al}"^ and by Le Hur and Seelig.^^ In this model, 
electron propagation between the two leads is excluded 
from the outset, accounting thereby for the strong sup- 
pression of clastic cotunneling at temperatures above the 
level spacing.^* The corresponding low-energy effective 
Hamiltonian is obtained from Eqs. (l)-(4) by setting 
iV = 1 and M = 2, and taking tLSm,i and 



t 



R 

Im 



tRSm,2- Converting for convenience from the 

channel labels to = 1 and m = 2 within the box to 
a = L and a = R, respectively, the resulting Hamil- 
tonian is given by 

a=L,R k,a 

a=L,R q,<y 

+ E E^«{4aCB,a.5-+H.C.}. (10) 
a=L,R k,q,a 

As in Eq. (9), the effective magnetic field h is equal to 
2Ec6N. 



1. Coherent propagation between the leads 

The Hamiltonian of Eq. (9) is equivalent in equilib- 
rium to the planner two-channel Kondo Hamiltonian in 
an applied magnetic field. This is best seen by first 
converting to a constant-energy-shell representation of 
the conduction-electron creation and annihilation oper- 
ators in each of the two leads and the quantum box, 
and then constructing generalized "bonding" and "anti- 
bonding" combinations of the two leads. Upon doing 
so the "anti-bonding" band decouples from the isospin 
S, while the "bonding" band plays the role of the single 
lead in Matveev's original mapping onto the planner two- 
channel Kondo Hamiltonian.^ At the degeneracy point, 
the capacitance of the single-electron transistor diverges 
logarithmically with decreasing temperature according to 
C oc (l/rR:)log(rRr/r), where Tk is the two-channel 
Kondo temperature:^^ 



ksTx = {Dy/gL + exp 



'iy/9L+9R. 



(11) 



Here D ~ 2Ec is the effective conduction-electron band- 
width. 



ga = Pa{0)PBiO)tl 



(12) 



{a = L,R) are the dimensionless tunneling conductances 
for the left and right junctions, and 



^T(e) = ^H£-'^jk) 



(13) 



are the underlying density of states (7 = L,R, or B). 
Throughout this paper we use the notation by which the 
argument of P'y{e) is measured relative to the correspond- 
ing chemical potential, i.e., Hl and iir for 7 = L and R, 
respectively, and /xb = for 7 = B. 

Away from equilibrium, it is still possible to construct 
"bonding" and "anti-bonding" combinations of the left 
and right leads, yet one can no longer dismiss the "anti- 
bonding" degrees of freedom as irrelevant. While the lat- 
ter degrees of freedom remain decoupled from the isospin 
S on the level of the Hamiltonian, they do couple to S 
through the effective density matrix, which assigns dif- 
ferent chemical potentials to the right and left leads. 
Hence both the "bonding" and "anti-bonding" combina- 
tions must be retained when computing the current for 
a finite bias. 



2. No coherent propagation between the leads 

Contrary to the Hamiltonian of Eq. (9), the Hamil- 
tonian of Eq. (10) corresponds in equilibrium to the 
planner four-channel Kondo Hamiltonian with channel 
anisotropy. Specifically, there are two distinct pairs of 
equivalent channels, one pair for each tunnel junction. A 
channcl-isotropic four-channel Kondo Hamiltonian is re- 
covered only for equal tunneling conductances for the left 
and right junctions, gL = gR- The latter conductances 
are given by 

5a =Pa(0)pSa(0)t' (W) 

(a = L,R), which differ from the expressions of Eq. (12) 
only in the separate density of states for the left and right 
modes within the quantum box: 



PBa(e) = ^^(e - CBka)- 



(15) 



For gi = gn = g, the low-energy physics of the tran- 
sistor is governed by the non-Fermi-liquid fixed point of 
the four-channel Kondo effect. The corresponding Kondo 
temperature is given by^^ 



ksTK = (Dg) exp 



(16) 



where D ~ 2Ec- Note that the exponent of Eq. (16) 
differs by a factor of \/2 from that of Eq. (11) under 
the same condition that 5l = 5fi = <?• Any asymmetry 
in the left and right tunneling conductances, gi, ^ gR, 
drives the system to a two-channel fixed point with one 
of the leads effectively decoupled from the box. 



III. NONCROSSING APPROXIMATION 

To obtain a reliable quantitative theory for the low- 
temperature transport at resonance for each of the mod- 
els of Eqs. (9) and (10), we resort to a recent adap- 
tation of the noncrossing approximation (NCA) to the 
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Kondo spin Hamiltonian with arbitrary spin-exchange 
and potential-scattering couphngs.^^ In this section, we 
formulate the Kondo-NCA (KNCA) for each of the mod- 
els of Eqs. (9) and (10), generalizing the KNCA to 
nonequilibrium. 

A. Slave- fermion representation 

To handle the isospin S, which we refer to hereafter as 
the impurity spin, we employ Abrikosov's slave-fermion 

representation.^"^ In this representation, one assigns a 
pseudo fermion to each impurity spin state according to 

/t |o) ^ = +1/2) , (17) 

/I |0) ^ \S, = -1/2). (18) 

This assignment corresponds to the replacement of the 
impurity spin operator by the bilinear pseudo-fermion 
operator 

where a are the Pauli matrices. The physical subspace 
corresponds to the constraint Nf = /j/^ = 1, which 
represents the fact that we are working within an en- 
larged Hilbert space. This constraint distinguishes the 
pseudo fermions from ordinary fermions. 

The advantage of the slave-fermion representation 
stems from the ability to use standard diagrammatic 
many-body techniques to calculate physical observables. 
The difficulty lies in implementing the constraint, which 
necessitates the introduction of a fictitious "chemical po- 
tential" A for the pseudo fermions. The latter is taken to 
minus infinity at the end of the calculation, as described 
below. For concreteness, let us focus in the following on 
the case where coherent propagation is allowed between 
the two junctions, i.e., the Hamiltonian of Eq. (9). The 
necessary modifications for the model of Eq. (10) are de- 
tailed in Appendix A. 

In conventional perturbation theory for a nonequilib- 
rium problem, one starts with an unperturbed system in 
equilibrium. All processes that drive the system out of 
equilibrium are then switched on adiabatically at some 
initial time to. For the problem at hand, one starts with 
two decoupled leads, each with its own chemical poten- 
tial. The unperturbed Hamiltonian Tio thus lacks the 
exchange term, and is given by 

^0 = X! + Ma)c],fe^Cafc<, 

a=L,R k,a 

+ E ^BA.a^B.a + ^ " ^)Afl- (20) 

g,a- 7=± 

Here e-y is equal to —7^/1, while A is a fictitious chemical 
potential for the pseudo particles. Note that Ho is bi- 
linear and diagonal in the single-particle operators cl^^^ 



and /|, which makes it a suitable starting point for di- 
agrammatic calculations. The initial density matrix of 
the system is also diagonal in the above single-particle 
operators, and has the form 

^° ^ Tr {e-/3(Wo-MiA'-L-i"KJVK)} ■ 

Due to the chemical potential A that was added to 
Ho, the statistical weight of the Nf = m subspace has 
an extra factor of e^^™ in Eq. (21). This allows one to 
project out the Nf = 1, physical subspace. Specifically, 
the average of any physical observable O can be expressed 
as 

(6)phys = ^ lim e-^^(diV/)A, (22) 

where 

Zi„,p= lim e-'^^{Nf)x. (23) 

A— ^ — oo 

Here subscripts A denote averages with respect to the en- 
larged Hilbert space. Since the Nf = subspace does not 
contribute to the averages of Eqs. (22) and (23) (other 
than through the normalization of po), then the leading- 
order terms in e^^^ come from the N f — 1, physical sub- 
space. The latter terms are the only ones to survive the 
A —00 limit. Note that, in practice, one can drop 
the Nf operator from the average of Eq. (22) for those 
physical operators O that give zero when acting on the 
Nf = subspace, which greatly simplifies the calcula- 
tions. We also note that .^imp corresponds in equilibrium 
to the "impurity contribution" to the partition function 
(see, e.g., Ref. 17). 

At to, the tunneling terms are switched on adiabati- 
cally, and the system evolves according to the full Hamil- 
tonian: Ho + Htun with 

ntun= E E^"{4aC^«<^/l/+ + H.C.}. (24) 

a—L.R. k.q.a 

After all transients have decayed, a new nonequilib- 
rium steady-state is reached, characterized by time- 
independent averages of physical observables such as the 
current operator. Within the enlarged Hilbert space, the 
steady-state average of such an operator O at time t = 
is given by 

(6)a= lim TT\poU^{0,to)6u{0,to)}, (25) 

to— > — oo L J 

where U is the timc-cvolution operator corresponding to 
the full Hamiltonian. Projection onto the physical sub- 
space is carried out according to Eq. (22). 
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FIG. 2. Diagrammatic representation of the KNCA 
approximation for the pseudo-fermion self-energies, for 
the Ifamiltonian of Eq. (9). Within the KNCA, the 
pseudo-fermion self-energies, Fig. (a), are approximated by 
a bubble with one bare conduction-electron propagator (full 
line) and one ladder propagator (dotted-dashcd line). The 
latter propagator is defined in Fig. (b). Full circles denote 
tunneling vertices. The ladder propagators are constructed 
from particle-hole bubbles, Fig. (c), with one bare conduc- 
tion-electron propagator and one fully dressed pseudo-fermion 
propagator (dashed line). There are two different bubbles: 
p1^', in which 7 = — and a = B, and p!^\ in which 7 = -|- 
and a is summed over L and R. The same convention for 
7 and a applies to the self-energy bubble of Fig. (a). Note 
that we include the matrix elements of the two edge vertices 
and a factor of ps(0)/p in our definition of p!^^\ while Pi^^ is 
defined without the matrix elements of the edge vertices, yet 
with an extra factor of 1/pb{0). Only particle-hole bubbles 
are included within the KNCA, while all other uoncrossing 
diagrams (e.g., the particle-particle bubbles) are omitted. 



B. Noncrossing approximation 

The key ingredients for the calculation of physical ob- 
servables are the pseudo-fermion Green functions. These 
include the retarded and advanced Green functions, 

t') = -i9{t - t') , (26) 

G<''\t,t')=ie{t'-t){{f,{t),f\{t')])^, (27) 

along with the lesser and greater Green functions, 

G<(^)(t,t') = (4(i').f7W)A, (28) 
G>''^\tX)^{m.f\{f))x. (29) 

Here curly brackets denote the anticommutator, and 
7 = ±. Once steady state is reached, the above Green 
functions regain timc-translational invariancc, and arc 
solely dependent on the time difference At = t — t' . It 
is therefore advantageous to switch over to the energy 



domain, by introducing the Fourier transforms with re- 
spect to At/h. In equilibrium, the lesser Green function 
is simply equal to the spectral part of the retarded Green 
function times 27r/(e), where /(e) is the Fermi-Dirac dis- 
tribution function. Away from equilibrium, when the 
effective distribution function is not known, both the re- 
tarded and lesser Green functions are explicitly needed 
in order to compute physical observables. 

In practice, the pseudo-fermion Green functions enter 
the calculation of physical observables in their projected 
forms, which read 



G!:(e)= lim G:(^)(e-A), 

A— — 00 

G^(e)= lim G"(^)(e-A), 
G<(e)= lim e-^^G<^\e-X), 

' A— >— 00 

G>(e)= lim G>(^)(e-A). 



(30) 

(31) 
(32) 
(33) 



These projections are analogous to the ones used in equi- 
librium, when G*^ coincides with the negative-frequency, 
or "defect", spectral function.^'' Note that, unlike G^ 
and G^, the lesser Green function G< has no contribu- 
tion from the Nf = subspace. Rather, its leading-order 
term comes from the Nf = 1 subspace, which decays to 
zero as e^^. The extra e~^'^ exponent in Eq. (32) is re- 
sponsible for canceling this decay to zero when the limit 
A — > — 00 is implemented. 

The projected Green functions have standard forms in 
terms of the projected self-energies. Specifically, Gl^'°'{e) 
and G^-'^(e) are equal to 



G!:'«(e) 



Sr.a / \ ) 
/ 7 \^) 



G<'>(e) = S<'>(e)|G!:(e)| 



(34) 
(35) 



where EI^'"(e) and 5]^'^(e) are obtained from their un- 
projected counterparts according to the projection rules 
of Eqs. (30)-(33). The KNCA^^ consists of a particular 
set of diagrams for the pseudo-fermion self-energies, de- 
picted in Fig. 2. Figure 2(c) shows the building block for 
the NCA self-energy diagrams of Fig. 2(a). It consists 
of a particle-hole bubble, with one fully dressed pseudo- 
fermion line and one bare conduction-electron line. The 
KNCA self-energy features a ladder of such bubble dia- 
grams. Fig. 2(b), with vertices inserted in between. Each 
ladder has an odd number of bubbles, as the incoming 
and outgoing pseudo-fermion lines share the same isospin 
label. For the Hamiltonian of Eq. (9), the physical spin 
is conserved along the ladder, whereas both the physical 
spin and the lead index (i.e., a = L,R) are conserved for 
the Hamiltonian of Eq. (10). 

Focusing on the Hamiltonian of Eq. (9), there are two 

separate ladders, labeled by the index of the outgo- 
ing and incoming pseudo-fermion lines. Projecting Dj^-* 
according to the rules laid out in Eqs. (30)-(33), one ob- 
tains 
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l-gPl{e)PL{ey 



gP<(e)+fl^P^(e)|P;(e) 
\l-gP^{e)PL{e)\^ 



(36) 



(37) 



where ql and qr are the dimensionless tunneUng conduc- 
tances for the left and right junction, defined in Eq. (12), 
g = gL + gR is the sum of the two dimensionless conduc- 
tances, and 



(38) 

(39) 
(40) 

(41) 



PU^) = - E T / G;(e')/(e'-e-/x„) 

X iya{e-e-iJ,a)de, 
P<{e) = -J2 - r G^(e')/(-e' + e+Ma) 

a=L,R ^ 

X Va{e'-€-fia)de\ 

/OO 
GL{e')f{e'-e)vB{e'-e)de', 
-OO 

/OO 
Gt{e')f{-e' + e)M^'-e)dt'. 
-OO 

Here v-yi^) — pj{e)/ pj{0) with 'y = L,R, or B is the 
reduced density of states. 

In Eqs. (36) and (37) we adopted the convention by 
which the -D_ ladder includes the matrix elements of the 
two end-point vertices, and is multiplied by pb(0). The 
_D+ ladder does not include; the matrix elements of the 
two end-point vertices, but is multiplied hy g/pB{0). We 
further restricted attention to the case of a zero physical 
magnetic field (not to be confused with h = 2EcSN), in 
which case all dependences on the physical spin index a 
drop. With these notations the KNCA self-energies take 
the form 



s;(e) 

S^(e) 
Sf(e) 



-2 E - r ^;(e')/(-e+e'+M„) 

X i/„(e-e'-^a)de', 

/OO 
D<(e')/(e-e'-M«) 

X Ua{€-e'-iJ,a)de', 

/OO 
DL{e')f{-e+e')Me-e')de', 
-OO 

/OO 
D<{e')f{e-e')Me-e')de', 
-OO 



(42) 

(43) 

(44) 



where an extra factor of two comes from summation over 
the two equivalent spin orientations. 

Equations (36)-(45) represent the complete summa- 
tion of the KNCA class of diagrams for the pseudo- 
fermion self-energies, in the case where coherent prop- 
agation is allowed between the two leads [i.e., the Hamil- 
tonian of Eq. (9)]. For gvB{e) = gLi^L{e) + gRVRie),"^^ 



these equations properly reduce in equilibrium to those 
of Ref. 22 for the planner two-channel Kondo Hamilto- 
nian.^*^ The adaptation of Eqs. (36)-(45) to the Hamil- 
tonian of Eq. (10) is specified in Appendix A. 



IV. FORMULATION OF THE CURRENT 

Our next goal is to formulate the current in terms of 
the pseudo-fermion Green functions introduced in the 
previous section. We begin our discussion with the model 
of Eq. (9). 



A. Coherent propagation between the two junctions 

The operator 7„, describing the electrical current flow- 
ing into lead a, is given by the time derivative of the 
charge operator for that lead, Qa- Within the slave- 
fermion representation of the Hamiltonian of Eq. (9) one 
obtains 

L = i^t^ E {4aCSfe'a/l/+ - H.C.} . (46) 

k,k' ,(T 

Applying the projection procedure of Eq. (22), the time- 
averaged current Ia{V) reads 



/a(V') = -2ji„^Im{e<,(t,i)}, 



(47) 



where Gaaitft') is equal to 



k,k' 



with 



Fka{t)=CBka{t)fUt)f+it)- 



(48) 
(49) 



In steady state, Q^^ {t, t') depends solely on the time 
difference At = t — t'. Switching over to the energy 
domain and applying standard perturbation theory one 
obtains the exact relation 

g<„{e) = t„ [GLp(e)r<(e-Ma) + Gf^^{e)K{e- n^)] , 

(50) 



(45) where GJ^ (e) is the Fourier transform of the response 



function 

GL^M = ^ lim e-^'Y.(^l'At')PkAt) 



k,k' 



(51) 



^imp(^) is its retarded counterpart, and r^(e) and rj^(e) 
are equal to 
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FIG. 3. The KNCA diagram for Gimp. Within the KNCA, 
Gimp reduces to a simple bubble, with one fully dressed /+ 
pseudo-fermion propagator (dashed line), and one £)+ ladder 
propagator (dotted-dashed line). 



Figure 3 shows the KNCA diagram for 



Sim- 



ilar to the standard NCA formulation of the impurity 
Green function for the Anderson impurity model, Gimp 
is given by a simple bubble diagram, which consists of 
one fully dressed /+ pseudo-fermion propagator and one 
ladder propagator. The resulting KNCA expressions 
for Gr^p(e) and Gf^^{t) read 



Pb(0) 



-G;(e + e')^^(e')]^e', (57) 



r<(e)=27rp„(e)/(e), 



y_oo e-e'-i 



Physically. T^a (e) 



Z(5 



(52) 



(53) 



is the conduction- 



electron T-matrix for scattering from lead a to lead a. 
Hence Gimp is analogous to the dressed impurity Green 
function in the Anderson impurity model. Indeed, insert- 
ing Eq. (50) into Eq. (47) yields 



?ipB(0) 



[2^Aimp(e).f (e - Ma) - G<,p(e)] 
•%.Va{t- lia)de, (54) 



which has the same structure as the expression for the 
current through an Anderson impurity. Here j4iinp(e) = 
— (l/7r)Im{Grjjp(e)} is the spectral part of G\^^. 

The energy integration in Eq. (54) is restricted in 
practice to an interval corresponding to the chemical- 
potential difference between the leads, broadened by the 
temperature. To the extent that one can neglect the en- 
ergy dependence of Pa(e) on that scale, it is possible to 
eliminate Gj^p from the expression for the current by us- 
ing the identity II = -Ir = {qrIl - 9LlR)/g, to obtain 



Il{V) 



with 



_ 2e r 



Teff(e)[/(e-ML)-/(e-Mii)]rfe (55) 



Teff(e) 



Pb(0)<7 ^ 



^imp 



(e) 



(56) 



Equation (55) is exact in the wide-band limit, when 
the energy dependence of pa{() can be neglected. It has 
the familiar form of an integral of an effective transmis- 
sion coefficient Tes(^) times the difference of two Fermi 
functions. Similar to the case of tunneling through an 
Anderson impurity, 7^ff(e) depends strongly on the tem- 
perature and bias at low energies. Our main task is to 
evaluate this effective transmission coefficient, which re- 
quires knowledge of Gl^^{e) and Gj^p(e). We calculate 
the latter functions using the NCA class of diagrams. 



G<.p(e) = r G<(e + e')Im{2);(e')}de', (58) 



where 



^imp= y ^[G^(^)+G'<(e)]^• 



(59) 



Hence G-mp(e) and Gi"^p(e) reduce within the KNCA to 
convolutions of the 7 = + pseudo-fermion Green function 
and the D+ ladder. 



B. No coherent propagation between the two 
junctions 

Formulation of the current for the Hamiltonian of 
Eq. (10) follows the same basic steps presented in the 
previous subsection for the Hamiltonian of Eq. (9). The 
sole difference in Eqs. (46)-(49) enters through the c^fev 
operators, which acquire an additional lead index. The 
latter index carries over to Fkc, and G^p{t,t'), whose 
definitions are modified according to 



{t) = CBkaa {t)f_{t)U{t), 



(60) 



^imp,a(^' ^ ) 



1™ ^-''T.{^lk'At')F»kAt))- 



k,k' 



(61) 



Contrary to Eq. (54) for the Hamiltonian of Eq. (9), each 
of the expressions for Il{V) and Ir{V) involve then a 



different Gimp,a function: 



Ia{V) = J ^[27rAmp.a(e)/(e-/Xa) - G5^p,Je)] 



X Vai^ — ^J■a)de. 



(62) 



Here ga are the dimensionless tunneling conductances of 
Eq. (14), and ^imp,c«(e) is the spectral part of G-mp.a- 

In general, Gi,np,L and Gimp,_R are two distinct func- 
tions, with no simple relation. As a result, one can no 
longer exploit current conservation to eliminate G^^^ 
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from the expression for the current, as was done in 
Eq. (55). Hence, there is no simple analog to the effec- 
tive transmission coefficient of Eq. (56) wlien no coherent 
propagation is allowed between the two leads. This is to 
be expected, as electrons do not truly propagate between 
the two leads in this case. 

As before, Gimp,L(e) and Gimp.nXe) are still given 
within the KNCA by the bubble diagram of Fig. 3, which 
features, however, a different ladder for a = L and 
a = R. Explicitly, one has 

GCnp,Je) = -^^^ r [G<{e')Dl4e' - e) 

/7r(7ct^imp J -oo 

-G\{e + e')D<^{e')\de', (63) 

/CXD 
G<(e + e')Im{D;+(e')}de'. 
-OO 

(64) 



imp J —OO 



In the following section, we present our numerical re- 
sults for the differential conductance as obtained from 
Eqs. (54) and (62) using the KNCA. 



V. RESULTS 

We have computed the differential conductance for 
each of the models of Eqs. (9) and (10), by first solv- 
ing the KNCA equations for the pseudo-fermion Creen 
functions and ladders, and then evaluating the current 
using Eqs. (54) and (62), respectively. The differential 
conductance G{T, V) = dl/dV was obtained by numeri- 
cally differentiating the current with respect to V, with- 
out resorting to the wide-band limit leading to Eq. (55). 
The most difficult step in the above procedure is solu- 
tion of the KNCA equations, which requires simultane- 
ous solution of the retarded and lesser pseudo-fermion 
and ladder functions. This was achieved by repeated 
iterations of the KNCA equations until convergence is 
reached. As is always the case with numerical solutions 
of NC A- type equations, the key to high- precision numer- 
ics lies in a well-designed grid of mesh points that scatter 
points more densely near the threshold and peaks of the 
relevant functions. To this end, we have used a combi- 
nation of linear and logarithmic grids. As a critical test 
for the precision of our numerical code, we have checked 
in all our runs that the pseudo-fermion spectral func- 
tions fulfilled the spectral sum rule to within one part 
in one thousand (0.1%). Extreme precision was required 
at very low temperatures and bias, when ksT and \eV\ 
were both much smaller than the corresponding Kondo 
temperature UbTk- 

Throughout this paper we focus our attention on the 
degeneracy point, ft, = 0, and set = — /^i? = 
We further use a semi-circular conduction-electron den- 
sity of states with half width D for the leads and the 



quantum box: Pa{i) = /'7(0)-\/l — {e/D)"^. Here 7 is 
equal to L,R or B for the Hamiltonian of Eq. (9), and 
L,R,BL or BR for the Hamiltonian of Eq. (10). Note 
that we restrict ourselves for simplicity to a single joint 
bandwidth D for the box and the leads. Although this 
need not be the case in real systems, the precise details of 
the band cut-offs should not affect the low-energy physics 
of interest here.^^ 



A. Coherent propagation between the two junctions 

We begin our discussion with the case where coherent 
propagation is allowed between the leads, i.e., the Hamil- 
tonian of Eq. (9). As discussed above, the latter Hamil- 
tonian reduces in equilibrium to the planner two-channel 
Kondo model with poJ± = Thus, the low-energy 

physics of the model is governed by the corresponding 
Kondo temperature Tk, which we extract by fitting the 
slope of the ln(T) diverging term in the isospin suscepti- 
bility 



Xisospin 



dh 



h=0 



to the Bethe ansatz expression'^^ 

1 



Xisospin (T') 



2QkBTK 



HTk/T). 



(65) 



(66) 



For g = + gji = 0.04, which corresponds to poJ± = 
0.4 in the planner two-channel Kondo representation, we 
obtain UbTk/D = 3.85 x IQ-^. 

As discussed at length in Ref. 22, the KNCA esti- 
mate of Tk deviates from the correct value of Eq. (11). 
However, as shown by comparison with the exact Bethe 
ansatz solution, ^^■^'^ the KNCA gives surprisingly accu- 
rate results for the temperature and field dependence of 
the magnetic susceptibility in the isotropic two-channel 
model, when T and HBg.jH are rescaled with the ex- 
tracted Kondo temperature. We expect a similar level of 
accuracy to hold for the nonequilibrium current in each 
of the models under consideration. 



1. Symmetric coupling 

Figure 4 shows the temperature dependence of the 
zero-bias conductance for different values oi gL = gn- As 
expected of Kondo-assisted tunneling, the conductance 
is enhanced with decreasing temperature according to a 
universal curve. Although Tk is varied by nearly three 
orders of magnitude, all curves collapse onto a single 
line. Slight deviations from universality are seen above 
T/Tk ~ 10 for the larger couplings, which may be due 
to the asymmetric line shape of the imaginary part of 
the conduction-electron T-matrix within the KNCA (see, 
e.g.. Fig. 12 of Ref. 22). The latter asymmetry is en- 
hanced at high energies as the coupling is increased. 
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T/Tk 

FIG. 4. Temperature dependence of the zero-bias con- 
ductance, for the model of Eq. (9) with different values of 
Ql = Qr = 5/2. The corresponding Kondo temperatures are 
equal to ksTK/D = 8.4 x 10"®, 1.25 x 10"*, 7.25 x 10"", and 
3.85 X 10"^ for g = 0.01, 0.017, 0.026, and 0.04, respectively. 
With decreasing temperature, the conductance rnoriotonically 
increases according a universal curve (full line). Below Tk, 
the conductance crosses over to a VT temperature depen- 
dence, in accordance with the expected power-law behavior 
of the two-channel Kondo effect. This is demonstrated in 
the inset, where a log-log plot of [G(0,0) - G{0,T)]/B ver- 
sus T/Tk is shown and compared to ^JT (Tk (dashed line). 
Here G(0, 0) and B were obtained separately for each curve 
from a fit to the ypT temperature dependence of Eq. (67). 
All values of G(0, 0) and B so obtained fall within the range 
G(0, 0)/i/e^ = 1-225 ± 0.005 and Bhje^ = 0.68 ± 0.015. The 
extracted values for G(0, 0)/i/e^ all fall within 1% from the 
exact zero-temperature KNCA value of 7r^/8. 

For over a decade of temperature around Tk, the en- 
hancement in G{0, T) is approximately logarithmic in 
temperature. Below Tk, the conductance crosses over 
to a a/T temperature dependence. As demonstrated in 
the inset of Fig. 4, the low-temperature conductance can 
be successfully fit to the form 

G{0,T) = G{Q,0)~B^ (67) 

with G(0,0)ft/e2 = 1.225 ± 0.005 and Bh/e^ = 0.68 ± 
0.015. Such a square-root temperature dependence of 
the conductance is expected of the two-channel Kondo 
effect. It stems from the square- root energy dependence 
of the imaginary part of the zero-temperature, zero-bias 
conduction-electron T-matrix,^'^^ or Aij^pie) in this con- 
text. A similar square-root temperature dependence of 
the conductance was also found for tunneling through 
a two-channel Anderson impurity,^'^ and measured for 
zero-bias anomalies in ballistic metallic point contacts. 
Note that the extracted zero-temperature conductance 
corresponds to a mostly open channel (61% of the full 



O 
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FIG. 5. The differential conductance G{V, T) = dl/dV ver- 
sus voltage bias, for the model of Eq. (9) with Ql = Or = 0.02 
and different temperatures. The corresponding Kondo tem- 
perature is equal to UbTk/D = 3.85 x 10~^. With decreasing 
temperature, a zero-bieis anomaly develops in G{V,T). For 
T < Tk, the anomaly acquires a \/V voltage dependence, 
which is rounded off for voltages on the scale of the temper- 
ature. For T = 0, a sharp cusp is left in G(V, 0) at zero bias. 

conductance), and falls within 1% from the exact zero- 
temperature KNCA value of G{Q,Q)h/e^ = -k"^ The 
latter figure is obtained by taking the zero-temperature 
limit of Eqs. (36)-(45) and (57) and performing a Miiller- 
Hartmann type of analysis, "^^ reproducing the results of 
Cox and Ruckenstein for the NCA treatment of the mul- 
tichannel Anderson model. 

It should be emphasized that the above enhancement 
of the zero-bias conductance for single-mode junctions 
differs markedly from the case of wide tunnel junc- 
tions considered by Schoeller and Schon.^ Omitting inter- 
mode mixing, these authors found that G'(0, T) van- 
ishes as 1/ ln(T) at low temperatures, in accordance with 
the vanishing- coupling fixed point of the infinite-channel 
Kondo Hamiltonian. However, as recently shown by 
Zarand et al.,^^ for realistic junctions with inter-mode 
mixing there exists an exponentially small crossover tem- 
perature T*, below which two-channel Kondo behavior 
prevails even for wide junctions. Hence our results should 
also apply to wide tunnel junctions at sufficiently low 
temperatures, provided electrons can propagate coher- 
ently between the two leads. 

Figure 5 shows the differential conductance for the 
symmetric coupling g^^ = gj^ = 0.02 and different tem- 
peratures. With decreasing temperature, a zero-bias 
anomaly develops in G{V,T), the height of which is plot- 
ted in Fig. 4. For T < Tk, the anomaly acquires a 
W voltage dependence, which is rounded off for volt- 
ages on the scale of the temperature. For T <C Tk, one 
is left with a sharp cusp in G{V,T) at zero bias. Sim- 
ilar to the VT temperature dependence of the conduc- 
tance, also the VV voltage dependence of the differential 
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° T/T^=0.01 
□ T/T^=0.03 
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(eV/k^T)^'' 

FIG. 6. The function H = [G(0,T) - G{V,T)]/AT^^^ ver- 
sus {eV/kBTy^^ , for the model of Eq. (9) with symmetric 
coupling. Here ql = gn = 0.02, while A is extracted from 
the leading square-root temperature dependence of the con- 
ductance: G(0, T) = G(0, 0) - AVT. For T < , all curves 
collapse onto a single scaling curve. Deviations from scaling 
are seen for T as low as O.STk- At higher temperatures the 
curves fan out, signaling departure from the scaling regime 
of the two-channel Kondo effect. The dashed line shows for 
comparison the corresponding scaling curve for a noninter- 
acting tunnel junction with a square-root energy-dependent 
transmission coefficient (see text). 



conductance stems from the anomalous square-root en- 
ergy dependence of the imaginary part of the conduction- 
electron T-matrix in the two-channel Kondo effect. 

Perhaps the most distinctive signature of the two- 
channel Kondo effect, though, lies in the scaling form 
of the differential conductance with eV/ksT . Following 
Ralph et al.,^"^ we plotted in Fig. 6 the scaling function 



H{V,T) 



G{0,T) ^ G{V,T) 

aVt 



(68) 



versus x = {eV/kBTy^^. where A is extracted from 
the leading temperature dependence of the conductance: 
G(0,T) = G(0,0) - AVT [i.e., A is related to B of 
Eq. (67) through A = B/^/rK\. For T <^ Tk, all 
curves collapse onto a single scaling curve, indicating that 
H{V,T) reduces to a function of the single scaling vari- 
able eV/ksT . The resulting scaling curve features two 
qualitatively different regimes: x ^ 1 and x ^ 1. For 
X < 1, the temperature cuts off the two-channel Kondo 
effect. Hence H{x) is proportional to x^, corresponding 
to a quadratic voltage dependence of the differential con- 
ductance. By contrast, for x » 1 the voltage bias serves 
as the cut-off energy, and H{x) depends linearly on x. 
The crossover between these two regimes takes place for 
X ~ 1. With increasing T there are deviations from scal- 
ing. These are significant for T as low as 0.3Tr-. At 
higher temperatures the curves fan out, signaling depar- 



ture from the scaling regime of the two-channel Kondo 
effect. 

Note that the scaling curve of Fig. 6 is quite similar to 
the one computed for tunneling through a two-channel 
Anderson impurity,^^ and the one measured for zero-bias 
anomalies in ballistic metallic point contacts. While not 
unexpected, this resemblance of scaling curves is by no 
means obvious, since the voltage bias couples differently 
to the impurity within the two-channel Anderson model 
of Hettler et al?^ and the Hamiltonian of Eq. (9). In 
Ref. 21, the spin- up and spin-down conduction-electron 
sectors undergo the same chemical-potential splitting, as 
the bias couples to the two spin sectors in an identi- 
cal manner. By contrast, only the lead electrons (the 
isospin-up sector) experience a chemical-potential split- 
ting within the Hamiltonian of Eq. (9), which breaks the 
equivalence of the two isospin sectors. As indicated by 
Fig. 6, this qualitative difference in the coupling to the 
voltage bias has only a moderate effect on the scaling 
curve, which is more shallow in Fig. 6 than in Ref. 21. 

That the shape of the scaling curve depends only mod- 
erately on the microscopic details of the system can be 
understood from comparison to a simple toy model, con- 
sisting of a noninteracting tunnel junction with a square- 
root energy-dependent transmission coefficient: T(e) = 
7q — 7i|e|^/^. In this case one can compute the scaling 
function of Eq. (68) exactly to obtain^^ 



with 



n£V/fcBT) 
F(0) ' 



dx. 



(69) 



(70) 



Note that Eqs. (69) and (70) are independent of both Tq 
and 7i. As seen in Fig. 6, the scaling curve of Eq. (69) has 
the same general structure as that of the Hamiltonian of 
Eq. (9), despite the fact that the noninteracting model 
does not account for the strong temperature and volt- 
age dependence of the effective transmission coefficient 
of Eq. (56). Thus, while the precise shape of the scaling 
curve does depend on the microscopic details of the sys- 
tem at hand, its general structure appears to be dictated 
by the square-root energy dependence of the imaginary 
part of the conduction-electron T-matrix at zero temper- 
ature and zero bias. 



2. Asymmetric coupling 

Thus far we have focused on symmetric coupling to the 
right and left leads. Below we consider the effect of an 
asymmetry in the coupling to the two leads. Specifically, 
we vary x = qr/ {ql + Qr) while keeping g = Ql + 9r 
fixed, as to maintain the same Kondo temperature Tk- 
Our results are summarized in Fig. 7. 
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FIG. 7. The differential conductance G{V, T) = dl/dV ver- 
sus voltage bias, for the model of Eq. (9) with different ratios 

of the couplings to the right and left leads, x — gn / (gR. + Ql)- 
Here g — 0.04 is kept fixed in all curves as to maintain the 
same Kondo temperature ksTx/D = 3.85 x 10~^. The tem- 
perature is equal to T/Tk = 0.1, while Gaym = l.Ole^ /h is the 
conductance for symmetric coupling at T/Tk = 0.1. In going 
from up to down, x takes the values x = 0.5,0.64,0.69,0.8, 
and 0.94. With increasing asymmetry, G{V, T) is reduced. 
Reduction of the zero-bias conductance is according to the 
standard formula 4a;(l— a;) [left inset: circles are the computed 
KNCA conductances, the dashed curve marks the parabola 
4x(l — x)]. As a function of bias, an asymmetry develops in 
the differential conductance. This cisymmetry is an artifact 
of the KNCA. It is emphasized in the right inset, where the 
scaling function [G(0,T) - G(y, T)]/^T^/^ with A extracted 
separately for each curve, is plotted versus cV/KbT. 

The most notable effect of asymmetric coupling is the 

reduction in the differential conductance with increas- 
ing asymmetry. As expected, the zero-bias conductance 
varies with x according to 4a;(l — x), which stems from 
the fact that, in equilibrium, Aimp(e) in Eq. (56) depends 
solely on the sum g = gL+9R, and not on the individual 
ga's. This effect is well captured by the KNCA (Fig. 7, 
left inset), which serves yet as another check for our nu- 
merical procedure. 

However, the KNCA produces an artificial asymmetry 
in the voltage dependence of the differential conductance. 
To see this we note that for a symmetric density of states, 
Pai^) = Pq(~£)i the particle-hole transformation c , 
{Sx, Sy, Sz) {—Sx, Sy, —Sz) convBrts V to —V in both 
the Hamiltonian of Eq. (9) and the initial density ma- 
trix of Eq. (21).^^ Since the current operator of Eq. (46) 
is antisymmetric under this transformation, one has the 
identity Ia{V) = —Ia(—V). Hence the differential con- 
ductance is a symmetric function of the bias, regardless 
of whether gi, and gn are equal or not. 

It is easy to verify that the KNCA equations preserve 
this symmetry of the differential conductance in the case 




FIG. 8. Schematic renormalization-group (RG) flow dia- 
gram for the Hamiltonian of Eq. (10) in equilibrium. De- 
pending on the ratio of ql to ga, the system flows to one of 
three fixed points. For symmetric coupling, ql = gR, the sys- 
tem flows to the four-channel fixed point denoted by A. For 
any asymmetry in the couplings, the system flows to one of 
two two-channel fixed points (B or C), at which the initially 
weaker ga vanishes, while the initially stronger coupling renor- 
malizes to the intermediate-coupling fixed-point value of the 
two-channel Kondo effect. Thus, the box is strongly coupled 
to one lead, but effectively decoupled from the other lead. 
The unstable weak-coupling fixed point is denoted by D. 

of symmetric coupling to the two leads. This stems 
from the fact that, within the KNCA, the voltage bias 
enters the effective transmission coefficient of Eq. (56) 
only through the E+ and P+ bubbles of Figs. 2(a) and 
2(c), respectively. Since the two leads enter E+ and 
P+ with equal weights for g/j ~ g^, and since inverting 
the sign of V simply corresponds to interchanging the 
roles of the two leads, then the resulting bubbles, and 
thus 7^ff(e), remain unchanged upon flipping the sign of 
V. This is no longer the case for asymmetric coupling, 
when the two leads enter each of E+ and P+ with dif- 
ferent weights. Indeed, as seen in Fig. 7, for gn ^ g^ 
the KNCA differential conductance acquires an asym- 
metry that increases with increasing x. Note that sim- 
ilar (though more pronounced) asymmetries in the dif- 
ferential conductance were obtained within the NCA for 
tunneling through a two-channel Anderson impurity,^^ 
where G{V, T) = G{—V, T) is not guaranteed by symme- 
try. 

B. No coherent propagation between the two 
junctions 

We proceed with the case where no coherent propaga- 
tion is allowed between the leads, i.e., the Hamiltonian 
of Eq. (10). For general couplings, this Hamiltonian re- 
duces in equilibrium to the planner four-channel Kondo 
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FIG. 9. Temperature dependence of the zero-bias con- 
ductance, G{V = 0,r), for the model of Eq. (10) with 
gi^ = gn = 0.017. Here ksTo/D = 4.7 x 10"^ is a charac- 
teristic energy scale of order the Kondo temperature. Similar 
to the two-channel case of Eq. (9), the conductance mono- 
tonically increases with decreasing temperature, crossing over 
from an approximate logarithmic temperature dependence at 
intermediate temperatures to a T^^^ behavior at low tempera- 
tures. The latter behavior is demonstrated in the inset, where 
a log-log plot of [G(0,0) - G(0,T)]/B versus T /To is shown 
(open circles) and compared to {T /Tq)^^^ (dashed line). Here 
G(0,0) = 0.547eV^ and B = 0.195e'^/h were extracted from 
a fit to the T^^^ temperature dependence of Eq. (71). 



model with channel anisotropy. Specifically, within this 
mapping one has poJ± = '^sfgZ for the two channels asso- 
ciated with the left junction, and p{\J^ = 2^Jlju. for the 
two channels associated with the right junction. Only 
for = gji is an isotropic planner four-channel Kondo 
model recovered, which shows the qualitative difference 
between symmetric and asymmetric couplings. For sym- 
metric coupling, the model flows in equilibrium to the 
four-channel Kondo fixed point, while for any asymme- 
try it flows to a two-channel Kondo fixed point where 
one lead — that with the weaker bare — is effectively 
decoupled. This renormalization-group (RG) flow dia- 
gram of the model is illustrated in Fig. 8. As emphasized 
by Furusaki and Matveev,"'^'^ the resulting conductance 
identically vanishes at T = for any gi^ ^ gn. 

From the discussion above it is clear that, unlike in 
the case of the Hamiltonian of Eq. (9) , there is no single 
energy scale that governs the low-temperature, low-bias 
transport properties of the system for all ratios oi gL to 
gii. If for symmetric coupling, = gii = g, transport is 
governed at low energies by the four-channel Kondo tem- 
perature of Eq. (16), for asymmetric coupling there is a 
new crossover scale for the onset of two-channel behavior. 
Below we explore both regimes. 
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FIG. 10. The differential conductance G{V,T) = dl/dV 
versus voltage bias, for the model of Eq. (10) with 
Ql = Qr. = 0.017 and different temperatures. Here kBTo/D 
is equal to 4.7 x 10~^. Similar to the case of the Hamilto- 
nian of Eq. (9), a zero-bias anomaly develops in G(V,T) with 
decreasing temperature. At low temperature, the anomaly 
acquires a V^^^ voltage dependence, in accordance with the 
onset of the four-channel Kondo effect. 



1. Symmetric coupling 

Figure 9 shows the temperature dependence of the 
zero-bias conductance, for gi, = g^ = 0.017. In the ab- 
sence of a precise procedure for extracting the KNCA 
Kondo temperature in the four-channel case, we use 
the scale UbTq = 4.7 x 1Q~^D as a rough estimate 
of Tk- This estimate builds upon the fact that the 
KNCA four-channel Kondo temperature is proportional 
to exp[— 2/po./_l] rather than exp[— 7r/2/9o which 
translates to j'^^'-'"*) qc exp[— 1/y^] in the notation 
of Eq. (16). We expect To and the actual KNCA Kondo 
temperature to be related by a factor of order unity. 

In accordance with the onset of the four-channel Kondo 
effect, the conductance is enhanced with decreasing tem- 
perature, approaching the zero-temperature value of 
G(0,0) = 0.547e2//i. Somewhat below To, the conduc- 
tance crosses over from an approximate logarithmic tem- 
perature dependence to a T^/^ power-law form. The lat- 
ter behavior is demonstrated in the inset to Fig. 9, where 
the low-T conductance is successfully fitted to 



G(0,T) = G(0,0)-S 



T 



1/3 



(71) 



with G(0,0) = 0.547eV/?, and B = Q.l^be^/h. 

The anomalous T^/^ temperature dependence of the 
conductance at low T is again a manifestation of the 
Ifl^/"* energy dependence of the imaginary part of the 
zero-temperature conduction-electron T-matrix in the 
four-channel Kondo effect. ^'^^ The same power law also 
characterizes the voltage dependence of the zero-bias 
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FIG. 11. The function H = [G(0,T) - G{V,T)]/ AT^^^ 
versus {eV/ksTY'^, for the model of Eq. (10) with 
Ql = Qr = 0.017. Here A is extracted from the 
leading T^/^ temperature dependence of the conductance: 
G(0,T) = G(0,0) - AT^'^. For T <C To, all curves converge 
to a single scaling curve. However, pronounced deviations 
from scaling persist down to much lower temperatures than 
for the Hamiltonian of Eq. (9) (compare with Fig. 6). The 
dashed line shows for comparison the corresponding scaling 
curve for a noninteracting tunnel junction with an \t\^^^ en- 
ergy-dependent transmission coefficient. 



anomaly in the differential conductance, which is plot- 
ted in Fig. 10. 

Similar to the model of Eq. (9), the differential con- 
ductance for the model of Eq. (10) also displays scaling 
with V/T, but with a different power law. Specifically, 
Fig. 11 shows a plot of the scaling function 



H{V,T) 



G{0,T) - G{V,T) 



(72) 



versus x = {eV/kBTy/^, where A is extracted from 
the leading T^/"^ temperature dependence of the conduc- 
tance: G(0,T) = G(0,0) - ArV3 [i.e., A is related to 

B of Eq. (71) through A = B/t]I\ For T <C Tq, all 

lines converge to a single scaling curve, indicating that 
-ff (V, T) indeed reduces to a function of the single scaling 
variable eV/kBT. However, the approach to scaling with 
decreasing T, as well as the fanning out of the curves with 
increasing T, is notably slower than for the Hamiltonian 
of Eq. (9) (compare with Fig. 6). The slower approach 
to scaling is consistent with the fact that the sub-leading 
energy dependence of the imaginary part of the T = 
conduction-electron T-matrix, which breaks scaling, is 
stronger for the four-channel Kondo effect^'^^ (of order 
|g|2/3^ The slower fanning out at higher T stems from the 
weaker T^/^ temperature dependence of the denominator 
in Eq. (72), as compared to Vt in Eq. (68). Thus, scal- 
ing with V/T is a sharper diagnostic for the two-channel 
Kondo effect that develops for the Hamiltonian of Eq. (9) 



than it is for the four-channel Kondo effect that develops 
for the Hamiltonian of Eq. (10) with gj^ = g^. 

As in the case of the Hamiltonian of Eq. (9), the general 
shape of the scaling curve of Eq. (72) is largely dictated 
by the |e|^/^ energy dependence of the imaginary part of 
the zero-temperature, zero-bias conduction-electron T- 
matrix, T„a(e) = t1,G\^^^^{e) . To see this, we plotted for 
comparison in Fig. 11 the corresponding scaling function 
for a simple toy model, consisting of a noninteracting tun- 
nel junction with an |e| ^/'^ energy-dependent transmission 
coefficient: T(e) = Tq — 7i|e|^/^. The exact scaling curve 
for this toy model remains given by Eq. (69), but with a 
slightly modified F function:^^ 



ny) = 



(73) 



Clearly, the noninteracting model lacks the complicated 
temperature and voltage dependences of Gfj^p^ and 
^imp a -'^l- (62). Nevertheless, the scaling curve of 
Eqs. (69) and (73) is quite similar to that of the Hamil- 
tonian of Eq. (10) with ql = Qr, as seen in Fig. 11. Hence 
the overall shape of the scaling curve for the Hamiltonian 
of Eq. (10) is largely dictated by the \e\^^^ energy depen- 
dence of j4imp,c« (e) at zero bias and zero temperature. 



2. Asymmetric coupling 

As discussed above, one expects the zcro-tcmpcratme 
conductancx; to vanish for any asymmetry in the cou- 
pling to the left and right leads. In the limit of both a 
large asymmetry and strong tunneling to one lead (i.e., 
a nearly open channel), this was indeed shown to be the 
case by Furusaki and Matveev.^^ For general couplings, 
one expects the system to initially flow towards the four- 
channel fixed point denoted by A in Fig. 8, before curving 
towards one of the two-channel fixed points, either B or 
C in Fig. 8. Depending on the degree of asymmetry, this 
flow may result in a nonmonotonic temperature and volt- 
age dependence of the differential conductance. Below we 
explore this scenario within the KNCA. 

Figure 12 shows the evolution of the differential con- 
ductance with increasing x — gn/igR + Ql)- Here we 
have fixed the coupling to the right lead at qr — 0.04, 
and varied ql from gL/gR = ^ down to gh/gR = 0.25. As 
a reference energy, we use the two-channel Kondo tem- 
perature corresponding to sole coupling to the right lead, 
which is the only relevant low-energy scale for gR^ gL- 
For gji = 0.04 and gL = 0, this two-channel Kondo tem- 
perature is equal to kBT2ch/D = 3.85 x 10~^. All curves 
presented were computed for T/T2ch = 0.01. An im- 
portant point to notice is that the KNCA respects the 
symmetry G{V,T) = G{-V,T) for the Hamiltonian of 
Eq. (10) with gL ^ gn, in contrast to the Hamiltonian of 
Eq. (9). This is evident from the symmetric differential- 
conductance curves of Fig. 12. 
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FIG. 12. The differential conductance G{V,T) = dl/dV 
versus voltage bias, for the model of Eq. (10) with differ- 
ent X = gn/igR + Ql). Here qr = 0.04 and T/Tsch = 0.01 
are fixed in all curves. As a reference energy, we use the 
two-channel Kondo temperature obtained for solo coupling to 
the right lead with qr = 0.04: kBT^A/D = 3.85 x 10"^ 
With increasing x, there is a reduction of the overall differen- 
tial-conductance signal, similar to the one seen in Fig. 7. A 
qualitatively different behavior is found near zero bias, whore 
the pealc for symmetric coupling {x = 0.5) first splits with 
increasing x, leaving only a broadened dip for x = 0.8. For 
X = 0.53, the initial enhancement of the differential conduc- 
tance with decreasing V is associated with an initial flow to- 
wards the four-channel fixed point (point A in Fig. 8). The 
dip at = signals an eventual flow towards the two-channel 
fixed point where the left is lead decoupled (point C in Fig. 8). 
At zero temperature, the conductance vanishes for all x > 0.5. 



The obvious effect of increasing x in Fig. 12 is to reduce 
the overall differential-conductance signal. This is to be 
expected, and is also seen in Fig. 7 for the Hamiltonian 
of Eq. (9). A more dramatic effect takes place near zero 
bias. Here the peak that forms for symmetric coupling 
(a; = 0.5) first splits with increasing x, leaving only a 
broadened dip for the larger values of x in Fig. 12. In 
particular, for x close to but larger than 0.5 (exemplified 
by a; = 0.53), a nonmonotonic voltage dependence of 
the differential conductance is seen near zero bias. This 
behavior is consistent with the notion of an initial flow 
towards the four-channel fixed point, before departing 
towards the two-channel fixed point with the loft lead 
decoupled. For larger asymmetries (i.e., larger x), the 
system no longer flows close to the four-channel fixed 
point, leaving only a monotonic dip in G{V,T). 

The emergence of a new low-energy scale for asymmet- 
ric coupling is clearly seen in Fig. 12. For symmetric cou- 
pling, X = 0.5, the four-channel Kondo temperature Ti^h 
is the only relevant low-energy scale. For small asymme- 
tries, represented by x = 0.53, a new low-energy scale 
Tdip <C T^ch emerges, corresponding to half the width 
of the zero-temperature dip that opens in G{V,T). The 
new scale Tjip steadily grows with increasing asymmetry, 
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FIG. 13. Temperature dependence of the diflerential con- 
ductance G{V,T) = dl/dV, for the model of Eq. (10) with 
gR/iQR + 9l) = 0.8. Here ql = 0.01 and gR = 0.04. 
kBT2ch/D = 3.85 X 10~^ is the two-channel Kondo tem- 
perature obtained for sole coupling to the right lead with 
qr — 0.04. With decreasing temperature, a sharp dip de- 
velops in the diflerential conductance a± V — 0. Specif- 
ically, G(0,T) is reduced by a factor of 25 in going from 
r/T2ch = 3 down to T/T2ch = 0.01, consistent with the flow 
to the two-channel fixed point where the left lead is decoupled 
from the box. Inset: an enlarged image of the dip at low bias. 



saturating at Tdip ~ T2ch for large asymmetries (see, e.g., 
X = 0.8 in Figs. 12 and 13). In this limit, no traces are 
left of the somewhat smaller four-channel Kondo temper- 
ature T4ch- Due to the small temperatures involved we 
are unable to systematically study the evolution of Tdip 
for small asymmetries. Nevertheless, by going to lower 
temperatures wc can estimate that Tdip for x = 0.58 is 
at least ten-fold smaller than for x = 0.8. 

Figure 12 was obtained for a finite temperature. At 
zero temperature, the conductance should vanish for all 
X ^ 0.5. Although we are unable to directly access the 
zero-temperature limit, we can confirm this scenario by 
studying the temperature dependence of the differential 
conductance. Figures 13 and 14 show our results for 
X = 0.8. For such a large asymmetry, we are able to 
go down to temperatures much smaller than Tdip. Upon 
going from T/Tzch = 10 down to T/Tsch = 5 x 10"'', a 
500-fold reduction is seen in the conductance. The lat- 
ter continues to drop steeply even for T/T2ch as low as 
0.001, in accord with a vanishing T = conductance. At 
low temperatures, G{Q,T) vanishes linearly with T, in 
agreement with previous predictions. ^^'^^ However, the 
linear-in-T regime is restricted to extremely low temper- 
atures. While Tdip is of the order of T2ch, the linear tem- 
perature dependence of the conductance only sets in at 
T/T2ch ^ 0.01, i.e., two orders of magnitude below Tjip. 
Hence the linear-in-T regime is unlikely to be accessible 
experimentally, for weak single-mode junctions. 
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FIG. 14. Temperature dependence of the eonductance, for 
the model of Eq. (10) with = 0.01 and gu = 0.04. The 
conductance monotonically decreases with decreasing temper- 
ature, approaching zero for T — > 0. Below T/T2ch ~ 0.01, the 
conductance vanishes linearly with T, and is well approxi- 
mated by G(0,T)/i/e^ « 0.77T/T2ch (dashed line). Note that 
the linear temperature dependence only sets in at a temper- 
ature much lower than the characteristic width of the dip in 
the differential conductance (compare with Fig. 13). 

VI. DISCUSSION 

In this paper, we studied the linear and nonlinear 
transport properties of a single-electron transistor at the 
degeneracy point. Focusing on weak single-mode junc- 
tions, two opposing scenarios were considered; one by 
which electrons can propagate coherently between the 
two leads, and the other whereby no electron propaga- 
tion is allowed between the leads. In the former sce- 
nario, both leads are assumed to couple to the same mode 
within the box. This amounts to a single-mode version of 
the model used by Konig et al.^ to describe wide tunnel 
junctions. In the second scenario, introduced by Funisaki 
and Matveev,^^ each lead is coupled to a separate mode 
within the box. At the degeneracy point, each of these 
models corresponds to a different planner multichannel 
Kondo Hamiltonian the two-channel Kondo Hamilto- 
nian in the case where electrons can propagate coherently 
between the leads, and the four-channel Kondo Hamilto- 
nian with channel anisotropy in the case where electron 
propagation is excluded. 

Generalizing the noncrossing approximation to these 
two particular nonequilibrium Kondo- type Hamiltonians, 
distinct signatures of the multichannel Kondo effect are 
seen in the differential conductance. Primarily, for sym- 
metric coupling, a zero-bias anomaly develops with de- 
creasing temperature, with a zero-temperature conduc- 
tance of order unity. Specifically, within the KNCA we 
obtain G(0,0) = 1.225eV/i and G(0,0) = 0.547eV/i, 
respectively, for the case where electrons either can or 
cannot propagate coherently between the leads. As a 
function of temperature, the conductance crosses over 



from a characteristic ln(T) temperature dependence at 
intermediate T to a power-law dependence: G(0, T) = 
G(0,0) — AT'^ . Here -q is the anomalous power law of 
the leading energy dependence of the imaginary part 
of the zero-temperature conduction-electron T-matrix in 
the appropriate multichannel Kondo effect: 77 = 1/2 for 
the two-channel Kondo effect, and = 1/3 for the four- 
channel Kondo effect. ^'^^ The same anomalous power law 
also characterizes the low-temperature differential con- 
ductance, which varies as G{V, T) = G(0, T) - DV^ for 
UbT <eV ^ keTK. 

Scaling of the differential conductance with eV/kBT 
is perhaps the clearest fingerprint of the multichannel 
Kondo effect that develops. ^^'^^ As shown in Figs. 6 and 
11, the function 

= 2(5^11^ (74) 

with fcsT, eV <C ksTK reduces to a function of the sin- 
gle scaling variable x = {eV /ksT)'^ . For \x\ < 1, H{x) 
is proportional to x"^^^, reflecting the quadratic voltage 
dependence of G(V, T) for eV < ksT. By contrast, H{x) 
is linear in |x| for |a;| > 1, which follows from the volt- 
age dependence of G{V,T) for eV > ksT. The crossover 
between these two markedly different regimes takes place 
for |x| ~ 1, i.e., when eV ~ ksT. 

As shown in the text, the shape of the resulting scal- 
ing curves can be largely understood within a simple 
toy model, consisting of a noninteracting tunnel junc- 
tion with an |e|'' energy-dependent transmission coeffi- 
cient: T(e) = Tf) — 7i|e|''. Although this simplified model 
lacks the the nontrivial voltage and temperature depen- 
dences of the conduction-electron T-matrix for the two 
Kondo problems of interest, it does reproduce the gen- 
eral structure of the two scaling curves. Hence the latter 
structures are largely dictated by the \e\^ energy depen- 
dence of the imaginary part of the conduction-electron 
T-matrix at zero temperature and zero bias. We note, 
however, that scaling with eV/ksT is a sharper diagnos- 
tic for the two-channel Kondo effect that develops in the 
presence of coherent electron propagation, as compared 
to the four-channel Kondo effect that develops for the 
Hamiltonian of Eq. (10) with symmetric coupling. 

The scaling curve of Fig. 6 is quite similar to the one 
obtained for tunneling through a two-channel Anderson 
impurity,^^ where a different coupling of the impurity to 
the voltage bias was considered. In the two-channel An- 
derson model of Hettler et al.,"^^ both the spin- up and 
the spin-down conduction electrons undergo the same 
chemical-potential splitting. By contrast, only the lead 
electrons (the isospin-up sector) experience a chemical- 
potential splitting within the Hamiltonian of Eq. (9), 
which breaks the equivalence of the two isospin sectors. 
Despite this qualitative difference between the two mod- 
els, the resulting scaling curves look qualitatively the 
same, reinforcing the moderate dependence of the shape 
of the scaling curve on the microscopic details of the 
model. 
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A clear distinction between the two scenarios under 
consideration is revealed when the coupling to the two 
leads is asymmetric. For the Hamiltonian of Eq. (9), the 
main effect of an asymmetry in the coupling is to reduce 
the zero-temperature conductance by the conventional 
factor of AgLgiil (gL+gn)"^- By contrast, for the Hamilto- 
nian of Eq. (10), the zero-temperature conductance van- 
ishes for any asymmetry in the couplings. As first ar- 
gued by Furusaki and Matvcev,^^ this unexpected result 
stems from the fact that the Hamiltonian of Eq. (10) flows 
for any asymmetry to a two-channel Kondo fixed point 
where one lead is decoupled from the box. Depending on 
the degree of asymmetry, one can obtain then nonmono- 
tonic differential-conductance curves, as seen in Fig. 12 
for X = 0.53. The nonmonotonicity reflects an initial 
flow towards the four-channel fixed point, before curving 
towards the eventual two-channel fixed point where one 
lead is decoupled (sec Fig. 8 for schematic RG trajec- 
tories). For large asymmetries, the system never flows 
close to the four-channel fixed point, leaving only a dip 
in the low-bias differential conductance. In accordance 
with previous predictions, ^^'^^ the zero-bias conductance 
vanishes linearly with T. However, as seen in Fig. 14, the 
linear-in-T regime is restricted to extremely low temper- 
atures, and is unlikely to be accessible experimentally for 
weak single-mode junctions. 

Throughout our discussion we have tactfully assumed 
that the level spacing inside the box is sufficiently small 
for a continuum-limit description to be used. Indeed, 
there are two prerequisites for observing a fully devel- 
oped Kondo effect: (i) the charging energy Ec must be 
sufficiently large for an experimentally accessible Kondo 
temperature to emerge, and (ii) the level spacing must be 
sufficiently small as not to cut off the multichannel Kondo 
effect. As emphasized by Zarand et al.,^'^ it is difficult to 
simultaneously fulfill these two conditions in present-day 
semiconductor devices (although some signatures of the 
two-channel Kondo effect were recently observed in the 
charging of a semiconductor quantum box^^). It is our 
hope that the host of signatures provided by this paper 
will assist in discerning the two scenarios proposed in the 
literature, once these experimental difficulties are over- 
come. 
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APPENDIX A: KNCA EQUATIONS IN THE 
ABSENCE OF COHERENT PROPAGATION 
BETWEEN THE LEADS 

Equations (36)-(45) were derived for the Hamiltonian 
of Eq. (9). In this appendix, we specify the corresponding 
KNCA equations for the Hamiltonian of Eq. (10). 

The main modification to the KNCA equations for the 
Hamiltonian of Eq. (10) comes from the ladder propaga- 
tors -Djf\ which carry an additional lead index a = L, R. 
Indeed, contrary to the Hamiltonian of Eq. (9), both the 
spin and lead indices arc conserved along each ladder for 
the Hamiltonian of Eq. (10). In the absence of an ap- 
plied magnetic field, when the two spin orientations are 
equivalent, the resulting ladders are independent of the 
spin index, but do depend on the lead index. 

The adaptation of Eqs. (36)-(41) to the Hamiltonian 
of Eq. (10) reads 



D: 



'a±{e) = — ^ 



-5a-P^+(e)Pg^_(e)' 



|l-5gP„V(e)Pg^_(e)|' 



(Al) 



(A2) 



where g^ and gu. are the dimensionless tunneling conduc- 
tances for the left and right junction, defined in Eq. (14), 
and 



PL+{^) 
P^+i^) 

P^e) 
PUe) 



/oo 
G;(e')/(e'-e-A^o) 
-OO 

X ^^a(e' — e — fia)de\ (A3) 

/oo 
G<(e')/(-e' + e + Ma) 
-OO 

X Ua{e' — e — /J-a)de\ (A4) 

/oo 
G''_{e')f{e' -e)uBaie' -e)de', (A5) 
-OO 

/OO 
G<{e')f{-e' + e)uBaie' -e)de'. (A6) 
-OO 



Here we have included the matrix elements of the two 

end-point vertices within the Da± ladders, and mul- 
tiplied Da+ and Da- by Pq(0) and pbq(O), respec- 
tively. Va{e) = Pa{e)/pa{0) and i^Ba(e) = PBaie)/ PBa{0) 
with a = L, R are the reduced density of states in the 
leads and in the box. The resulting pseudo-fermion self- 
energies acquire the form 



E;(e) = -2 ^ r I?;+(e')/(-e + e' + m«) 

a=L,R-^-°° 

X Uaie - e' - fj,a)de' , (A7) 

/OO 
D<^{e')f{e-e'-na) 
_ -OO 



a=L,R'' 



17 



y. i'a{e - e' - fict)de' , (A8) 
S!l(e) = -2 ^ r Dl_{e')f{-e+e') 

X UBa{e - e')de' , (A9) 
S<(e) = -2 J2 r D<_{e')f{e-e') 



a=L,R 



X i^Ba(e — e')de' 



(AlO) 



where the extra factor of two comes from summation over 
the two equivalent spin orientations. 

Equations (Al)-(AIO) are the complete summation 
of the KNCA class of diagrams for the pseudo-fermion 
sclf-cncrgics, in the case where no coherent propagation 
is allowed between the leads [i.e., the Hamiltonian of 
Eq. (10)]. They are analogous to Eqs. (36)-(45) for the 
Hamiltonian of Eq. (9). 
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